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ABSTRACT 

Supermassive black holes in giant elliptical galaxies are remarkably faint given their expected accretion 
rates. This motivates models of radiatively inefficient accretion, due to cither ion-electron thermal 
decoupling, generation of outflows that inhibit accretion, or settling of gas to a gravitationally unstable 
disk that forms stars in preference to feeding the black hole. The latter model predicts the presence of 
cold molecular gas in a thin disk around the black hole. Here we report Submillimeter Array observations 
of the nucleus of the giant elliptical galaxy M87 that probe 230 GHz continuum and CO(J=2-l) line 
emission. Continuum emission is detected from the nucleus and several knots in the jet, including 
one that has been undergoing flaring behavior. We estimate a conservative upper limit on the mass 
of molecular gas within ~ 100 pc and ±400 kms -1 line of sight velocity of the central black hole of 
~ 8 x 1O 6 M , which includes an allowance for possible systematic errors associated with subtraction of 
the continuum. Ignoring such errors, we have a 3<r sensitivity to about 3x 10 6 M Q . In fact, the continuum- 
subtracted spectrum shows weak emission features extending up to 4a above the RMS dispersion of the 
line- free channels. These may be artifacts of the continuum subtraction process. Alternatively, if they 
are interpreted as CO emission, then the implied molecular gas mass is ~ 5 x 10 6 M Q spread out over a 
velocity range of 700 kms -1 . These constraints on molecular gas mass are close to the predictions of the 
model of self-gravitating, star- forming accretion disks fed by Bondi accretion (Tan & Blackman 2005). 

Subject headings: accretion — stars: formation — galaxies: active, individual (M87), jets 

1. INTRODUCTION 

The most massive black holes in the Universe reside in the centers of giant elliptical galaxies. The (3.4 ±1.0) x 10 9 M Q 
(Macchetto et al. 1997) behemoth in M87 in the Virgo Cluster is one of the nearest (d = 16 ± 1.2 Mpc, Tonry et al. 
2001, for which 1"=78 pc). Such black holes are thought to have accumulated most of their mass by accretion of gas from 
thin disks, which efficiently radiate about 10% of the rest mass energy of the gas. These systems are likely to power the 
luminous quasars that are common at high redshifts. 

Most black holes we observe today in the local Universe are not surrounded by very luminous accretion disks. This is 
surprising because these black holes are often embedded in gas that should be able to accrete. In the centers of elliptical 
galaxies the gas is typically at temperatures of kT ~ 1 keV, i.e. T ~ 10 7 K, corresponding to sound speeds of 400 kms -1 , 
and at particle densities of about 1 cm -3 . These properties can be measured relatively accurately with high-spatial- 
resolution X-ray telescopes. For example the Chandra X-ray Observatory measures kT = 0.8 keV and n e — 0.17 cm -3 in 
the central kpc of M87 (Di Matteo et al. 2003) . The expected Bondi accretion rate to the black hole is given approximately 
by the product of the gas density, sound speed, and area of a sphere with a radius such that the escape speed at this 
distance is the sound speed. For M87 the expected accretion rate is about 0.04 M Q yr _1 and the luminosity of a thin 
accretion disk would then be 2 x 10 44 ergs s -1 (Tan & Blackman 2005, hereafter TB05). However, the radiative luminosity 
of the nucleus is observed to be only 10 42 ergs s _1 (Birctta, Stern & Harris 1991). 

It is possible that the radiative luminosity is only a fraction of the total accretion luminosity since accretion likely 
powers the observed M87 jet. The mechanical luminosity of the observed M87 jet has been estimated to be as large as 
2 x 10 43 ergss -1 (Reynolds et al. 1996) by assuming a minimum pressure in the 5 kpc radio lobes to obtain the total 
energy, in combination with synchrotron spectral aging to obtain a lifetime of 10 6 years. These estimates are sensitive to 
any uncertainties in the aging. Also, it is not certain that the M87 outflow luminosity has been steady from the beginning 
of the radio lobe formation to the present. Nevertheless, if the jet luminosity were steady and the above estimate correct, 
then the discrepancy between a model involving the extraction of 10% of the rest mass energy of an accretion flow equal 
to the present Bondi Accretion rate and the total output from accretion (radiative + mechanical) would still be a factor 
of 10. 

To help explain these discrepancies, models of radiatively inefficient accretion flows (RIAFs) have been proposed. In 
the simplest variant, the cooling time of the gas is long compared to the accretion time, so much of the energy is advected 
into the black hole (Ichimaru 1977; Narayan & Yi 1995; Quataert & Gruzinov 1999). In other variants, the inability 
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of the gas to cool leads to reduced accretion rates near the black hole due to the generation of outflows (Blandford & 
Begelman 1999) or convective motions (Quataert & Gruzinov 2000). A qualitatively different model allows the gas to cool 
to form a thin disk, but prevents accretion of the gas by having it form stars (TB05). These stars are then on relatively 
dissipationless orbits about the black hole. A prediction of this model is the presence of a star-forming disk inside the 
Bondi radius of the black hole. Since all known star formation in the Universe occurs in cold, molecular gas, that is best 
traced by the rotational transitions of carbon monoxide (CO), this model predicts that CO emission should be seen from 
such a circumnuclear disk. 

In the case of M87, TB05 predicted a molecular gas mass of about (1 — 5) x 10 6 M Q , assuming the disk is self-regulated 
to have a Toomre Q parameter of unity by the energy input from star formation at a rate equal to the Bondi accretion 
rate. This star formation rate is consistent with the observed Ha luminosity of the nuclear disk, although it is possible 
that the AGN also contributes significantly to the ionization of this gas. The quoted variation in this estimate is due to 
whether the disk is optically thick or thin to the heating radiation from the stars. Additional systematic uncertainties at 
the factor of a few level are caused by uncertainties in the stellar initial mass function (IMF), radial distribution of the 
gas and star formation, and the black hole mass. 

The deep potential well of the M87 black hole causes large orbital speeds in the disk. Given the above mass estimate, 
these are equal to 382 ± 60kms~ 1 at a distance of 100 pc, about the extent of the observed Ha (Ford et al. 1994) and 

II disk (Macchetto et al. 1997). From the latter study, the best fit value of the inclination angle of the disk's rotation 
axis to our line of sight is 51°, the systemic velocity is 1290 kms -1 , and the observed FWHM along a slit running through 
the center of the disk and 1" in extent is about 700 km s _1 . 

If CO emission emerges from the same region, then one expects a similarly broad line, which presents a problem for its 
detection since most radio telescopes have relatively narrow band passes. The Submillimeter Array (SM^Q, Ho, Moran, 
& Lo 2004) has a large band pass of about 2 GHz, corresponding to 2600 kms -1 at the frequency of CO(J=2-l) line 
emission (i.e. 230 GHz, 1.3mm). Having significant line-free regions in the bandpass is important for determining the 
continuum level, since the immediate vicinity of the black hole is known to be a strong (~1.5 Jy) radio continuum source 
at A=l mm. For these reasons, we used the SMA to perform a search for molecular gas in the nuclear disk of M87. 

1.1. Previous searches for CO in M87 

Braine & Wiklind (1993) searched for CO in M87 with the IRAM 30m single-dish telescope, with a beam size of 
13"(=1.0 kpc) for the CO(2-l) search at 230 GHz and 21" (=1.6 kpc) for the CO(1-0) search at 115 GHz. Thus these 
observations did not resolve the Ha circumnuclear disk. Another limitation was the relatively narrow bandpass, cor- 
responding to 700 km s -1 for CO(2-l) and 1300 km s -1 for CO(1-0). Braine & Wiklind (1993) derived an upper limit 
to the molecular gas mass of 3 x 10 6 Mq from their CO(2-l) data (based on a 3cr uncertainty of 0.9 Kkms -1 over 
the assumed 500 km s _1 line width, a J=2-l to J=l-0 line ratio of unity, a standard CO to H2 conversion factor of 
N (R2) / Ico(i-o) — 2.3 x 10 20 mol cm _2 /(K kms -1 ) for optically thick emission at temperatures ~ 10 — 20 K, and a 
geometric area of ~ 7r(500 pc) 2 ) and a limit of 1 x 10 7 Mq from their CO(1-0) data. The limit from CO(2-l) is not valid 
for a line that is broader than about 500 kms -1 . Braine & Wiklind (1993) also derived more stringent limits assuming 
optically thin CO emission, although this would require iV(H2) < 10 21 cm -2 for typical molecular cloud conditions, which 
would not be expected if the gas was organized into self-gravitating molecular clouds, as in the model of TB05. 

Combes, Young & Bureau (2007) also searched for CO(1-0) and CO(2-l) emission from the central region of M87 with 
the IRAM 30m single-dish telescope. From the lack of detected CO(1-0) emission in their 23"(1.8 kpc) FWHM beam and 
1300 kms -1 bandpass, they estimate a 3er upper limit of 1.1 x 10 7 Mq of molecular gas (assuming N / Ico{i-a) 
'.' x 10 20 mol cm~ 2 /(K kms -1 ), a 300 kms" 1 line width, and a geometric area of ~ 7r(900pc) 2 ). Using similar observations 
of CO(1-0) with the IRAM 30m telescope, Salome & Combes (2008) reported an upper mass limit of 7.7 x 10 6 M Q (la 
corrected to our adopted distance of 16 Mpc) evaluated assuming a 300 kms -1 line width. We will typically consider 3a 
upper limits, which for their results is 2.2 x 1O 7 M . 

2. OBSERVATIONS 

The M87 nucleus was observed with the SMA at 230 GHz in the extended configuration with 7 antennas in the array 
during the night to February 10th 2006. The primary beam at the given frequency was ~ 55", and the unprojected 
baselines ranged between 28 and 226 m. The phase center of the observations was R.A. (J2000.0) 12™30 m 49."4 and Dec. 
(J2000.0) 12°23'28."0. We had excellent weather conditions with measured zenith opacities r(225GHz) between 0.03 and 
0.06. 

Bandpass calibration was done with the quasar 3C273. Based on the spectrum taken toward the additionally observed 
quasar 3C279, we estimate the bandpass accuracy to be good to approximately 1%, which was only possible because of 
the excellent weather conditions. The flux calibration was performed with measurements of Callisto. We estimate the 
absolute flux uncertainties to be within 20%. Phase and amplitude calibration were done via frequent observations of the 
quasar 3C273, about 10.3° from the phase center. The absolute positional accuracy is estimated to be good to within 0.1". 
The receiver operated in a double-sideband mode with an IF band of 4-6 GHz so that the upper and lower sidebands were 
separated by 10 GHz. The central frequency of the upper sideband was centered at the CO(2-l) line at 230.538 GHz, the 

1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 
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v\ sr was set to 1200 km s -1 . The correlator had a bandwidth of 2 GHz and the channel spacing was 0.8125 MHz. Measured 
double-sideband system temperatures corrected to the top of the atmosphere were between 70 and 200 K, depending on 
the zenith opacity and the elevation of the source. The initial flagging and calibration was done with the IDL superset 
MIR originally developed for the Owens Valley Radio Observatory (Scoville et al. 1993) and adapted for the SMA3- The 
imaging and data analysis were conducted in MIRIAD (Sault, Teuben, & Wright 1995). 

Since the upper limit of the CO(2-l) line intensity is < 5% of the mm continuum, we averaged the whole upper 
sideband bandpass to produce our mm continuum image. This small line-to-continuum ratio made the CO line imaging 
and the connected continuum subtraction to the data a very difficult task. We tried various different approaches for the 
continuum subtraction using the MIRIAD and the AIPS software packages to verify consistent results. The final dataset 
was produced in MIRIAD using the task UVLIN and setting a spectral window between 700 and 1700 km s" 1 to ensure 
exclusion of any potential line contamination in the continuum subtraction fits. The final synthesized beams for the 
continuum and line datasets are 1.2" x 0.8". Our sensitivity was dynamic-range limited by the side-lobes of the strongest 
emission peaks. The la rms dispersion of the continuum data in regions of the spectrum expected to be free of CO(2-l) 
line emission and averaged over 100 kms -1 channels is 6.8 mJy beam -1 . Since our shortest baseline was 28m, this dataset 
filters out all emission on scales larger than ~ 12". 

3. RESULTS 

3.1. 230 GHz Continuum Image of Nucleus and Jet 

The 230 GHz continuum image of the M87 nuclear region is shown in Fig. [TJ together with the Ha image from HST 
(Ford et al. 1994). At 230 GHz, the nucleus is detected at a position (0.3", 0.0") relative to the phase center (SgJ). The 
peak and integrated fluxes of the nucleus are 1520 mJy beam -1 and 1770 mJy, respectively. The HST image was shifted 
by 1.06" (consistent with HST pointing uncertainties) to the south-west to align the nucleus and jet axis with the sub-mm 
data. 

Figure [1] also shows extended emission from the jet at about 1-2" (Knots HST-1 and HST-2) (Biretta, Sparks, & 
Macchetto 1999) and at 11-18" (Knots A, B, C) away from the nucleus (see also Fig. 1 of Perlman et al. 2001). Figure^] 
shows a profile of this emission along the jet axis and Table Q] lists the knot fluxes. Note that knots A and B are really one 

2 The MIR cookbook by Charlie Qi can be found at |http:/ /cfa-www. harvard.edu/~cqi/mircook. html 




Fig. 1. — 230 GHz continuum image of the nucleus and jet of M87 (contours), superposed with the optical HST image from Ford et al. 
(1994). The contouring of the 230 GHz data is from 18 to 90mjy beam -1 in 18mjy beam" 1 steps, then each contour level is at double the 
intensity of the previous. White and black contours are used for display purposes and are part of the same sequence. 
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physically-connected structure, and that knot B has two separate sub-peaks. When deriving fluxes for knots A and B we 
define their boundary at the location of the minimum flux between the knots, at about 12.5" from the nuclear emission 
peak. We assign 20% uncertainties to the knot fluxes (32). 

We compare our nuclear and knot fluxes to those measured in the infrared at 10.8 /im (2.78 x 10 13 Hz) by Perlman et 
al. (2001, their Table 1) and those measured at 3.4 mm (88.26 GHz) by Dcspringre et al. (1996, their §3.2, assuming 20% 
flux uncertainties): assuming a power law spectrum, F v oc v~ a , we derive spectral indices ai r _ mm and a mm , respectively 
(Table [T]). Care must be taken in the interpretation of these spectral indices since the nucleus and jet components may 
be variable: there is a 4.8 year time baseline between our 1.3 mm observations and the 10.8 /im observations and about 
a 12.5 year time baseline for our comparison with the 3.4 mm observations. For example, Steppe et al. (1988) found the 
nuclues of M87 is variable by about 1 Jy yr _1 at 90 GHz. Braine & Wiklind (1993) measured a 1.3 mm nuclear (< 13") 
flux of 1.0 Jy from their observations taken in Nov. 1991, a factor of 0.56 smaller than the value measured in this paper. 

3.1.1. Nucleus: Thermal Emission? 

Our value of ai r _ mm ~ 0.97 ± 0.04 for the nucleus is similar to the value of 0.88 predicted by the continuous injection 
synchrotron model (Heavens & Meisenheimer 1987) that was fitted to radio, optical and x-ray data by Sparks, Biretta, 
& Macchetto (1996) and Marshall et al. (2002). The spectral index we derive from comparison to the 3.4 mm data of 
Despringre et al. (1996) is a mm — 0.06 ± 0.26. This is similar to the index of 0.10 ± 0.02 derived from radio (1.47, 4.89, 
15.0 GHz) data by Biretta et al. (1991), but shallower than the index derived by comparing the 3.4 mm data to the radio 
data. 

We show the continuum spectrum of the nucleus in Fig. [3J The star-forming accretion disk model of TB05 predicts 
that there should be thermal emission from dust in the disk and the simple fiducial models of this emission are also shown 
in Fig. [3J These models assumed full reprocessing of the luminosity produced by stars forming (with Salpeter initial 
mass function down to 0.1 M©) at rates equal to the range of predicted Bondi accretion rates (0.036-0.15 M© yr -1 ). 
The luminosity was assumed to be radiated with a single temperature blackbody spectrum from a disk with radius of 
100 pc. In reality, only a fraction of any stellar luminosity is likely to be reprocessed by dust and there will be a range 
of temperatures extending to hotter values. Thus the flux in the Raylcigh- Jeans part of the spectrum (e.g. at 230 GHz) 
predicted by these models will tend to be an overestimate. More accurate models of the thermal emission are in the 
process of being developed (Tan, in prep.). 

In any case, our observed flux at 230 GHz is greater than the amount predicted by the TB05 models. The relatively 
shallow spectral index that we find from a comparison with the 3.4 mm data of Despringre et al. (1996) could be explained 
by the presence of a star- formation induced thermal component. Alternatively, it could be due to source variability or 
additional uncertainties in the measured fluxes. Note, Despringre et al. (1996) found extended (12.7" by 6.5") 3.4 mm 
emission around the nucleus with a flux of 1.4 Jy. Although this emission is aligned in a direction perpendicular to the jet, 
its scale is somewhat larger than that expected by a Bondi-fed accretion disk. Our 230 GHz interferometric observations 
do not show any evidence for such an extended component, although they lose sensitivity for scales > 12". 

Perlman et al. (2007) have found evidence for a thermal dust component in the nucleus of M87 based on Spitzer 
infrared data from 5 to 70 /Ltm. They fit models with temperatures in the range 45 to 65 K (for which the luminosities 
are 8.1 x 10 5 — 7.5 x 10 4 L© and the disk masses are 2.2 x 10 4 — 2.3 x 10 3 M©, respectively). Their data do not exclude 
the possible presence of additional cooler dust components. 

3.1.2. Jet Knots 

Our values of ai r _ mm ~ 0.6 for the jet knots A and B are very similar to the value of 0.68 predicted by the synchrotron 
models of Jaffe & Perola (1973), Kardashev (1962) and Pacholczyk (1970) presented by Marshall et al. (2002) and 
Perlman et al. (2001) based on the radio and optical data presented by Sparks, Biretta, & Macchetto (1996). Our value 
of ai r _ mm ~ 0.36 ± 0.05 for knot C is somewhat shallower than predicted, but it is the weakest of the knot components 
that we detect and we may have underestimated the uncertainty in its flux. From the comparison with the 3.4 mm data, 
we find a similar spectral index as the synchrotron models for knot A, ~ 0.6, but much steeper indices > 1 — 2 for knots B 
and C. The knot components were not well resolved by Despringre et al. (1996) so we evaluate the spectral index derived 
from the total flux of A, B and C: in this case we find a mm — 1.2 ±0.26, i.e. still significantly steeper than the synchrotron 



Table 1 

Nucleus and Jet Knot Fluxes 



Component 


1.3 mm Flux (mjy) 


C^ir— mm 




Nucleus 


1770 ± 350 


0.97 ±0.04 


0.06 ±0.26 


HST1/2 


180 ± 36 


n/a 


n/a 


A 


101 ± 20 


0.58 ±0.04 


0.55 ±0.26 


B 


98 ±20 


0.57 ±0.04 


1.28 ±0.26 


C 


23 ±5 


0.36 ±0.05 


2.4 ±0.26 
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Fig. 2. — 230 GHz continuum emission, integrated across a 4"wide strip, as a function of angular distance from the nucleus along the jet 
axis. 
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Fig. 3. — Continuum spectrum of the M87 nucleus. Data points are (from low to high frequency): 1.5, 5, 15 GHz (Biretta ct al. 1991); 
89 GHz (Despringre et al. 1996); 230 GHz (this paper); IRAS 100 /im, Spitzer-MIPS 70 and 24 /im and Spitzer-IRS (5-35 /im, dotted 
line) (Perlman et al. 2007); Gemini 10 /im (Perlman et al. 2001); optical and UV data (Sparks et al. 1996). The dashed and dotted lines 
show models of thermal dust emission from TB05, assuming luminosities of 5.9, 1.6 X 10 s Lq and temperatures of 15.5, 11.1 K, respectively. 
Note, these are simple, single-temperature, blackbody models that assume the entire luminosity due to star formation at rates equal to 0.15, 
0.036 Mq yr _1 with Salpeter initial mass function down to 0.1 Mq is reprocessed by dust in a 100 pc radius disk. 



models. This discrepancy could be due to variability in the knot fluxes or additional systematic uncertainties in the flux 
measurements . 

The peak and integrated fluxes of the inner sub- mm knot (sum of HST-1, known from optical and UV observations to 
be 0.85" from the nucleus, and HST-2) are 101 mJy beam -1 and 180 mJy, respectively. Knot HST-1 has been undergoing 
a flare in the last several years: its 2 keV x-ray flux increased by about a factor of 50 (Harris et al. 2006), peaking at a 
date of about April 2005 and declining to about half its peak value by August 2005. The rising part of the flare is also 
seen in UV (220 nm) and radio (15 GHz) data (Harris et al. 2006). It is likely that our measurement of the 230 GHz 
continuum in February 2006 includes a contribution from the declining part of this flare. 

3.2. Search for CO (J =2-1) Emission 

The continuum-subtracted spectrum extracted at the position of the continuum peak is shown in Fig. [4] Weak (~ 3 — 4er) 
emission features are seen in individual channels at 950 and 1350 kms , with 21 and 26 mJy, respectively. Note the 
dispersion of the line-free regions of the spectrum is 6.8 mJy in a 100 km s -1 channel. 

The continuum level is at 1.5 Jy, so the uncertainties in line detection are likely to be dominated by the systematic 
uncertainties associated with the continuum subtraction. The fact that the emission features are not in adjacent velocity 
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channels tends to argue in favor of their origin from systematic errors in the continuum subtraction process, although if 
the gas were distributed evenly in an annulus from 50 to 100 pc from the black hole, a double-peaked profile separated 
by about 500 kms -1 would be expected (Fig. [5]). If the molecular gas is further organized into spiral arms or a small 
number of discrete clouds, this could also lead to separated emission peaks. 

3.2.1. Applicability of Standard CO — H2 Conversion Factors 

The intensity of the CO line depends on the mass of molecular gas in the disk and its excitation conditions. If the 
gas is optically thick in the CO J=l-0 line then the column of gas is related to the intensity of the line via iV(H2) = 
1.8 ± 0.3 x 10 20 (/co,J=i-o/K kms _1 )cnr 2 (Dame, Hartmann, & Thaddeus 2001), based on observations of molecular 
clouds in the Milky Way. 

The applicability of this standard conversion factor to the nucleus of M87 is uncertain. There is some evidence that the 
conversion factor varies in different galaxies: e.g. in the center of M51, Israel , Tilanus & Baas (2006) find the CO-to-H2 
ratio is about a factor of 4 lower than the above Galactic value, i.e. four times less H2 for a given CO line intensity. 

X-ray ionization from the nucleus may be an important factor in determining CO abundance (Maloney, Hollcnbach, & 
Tielens 1996). The structure of these X-ray dissociation regions is controlled by the ratio of the x-ray flux to the gas density, 
described in the models of Maloney et al. (1996) by the dimensionless effective (attenuated) X-ray ionization parameter, 
£ c ff = 1-26 x 10~ 2 FxH^ 1 A^ ' 9 , where F x is the normally incident X-ray flux in ergscm _2 s _1 , n = 1000n3cm -3 is the num- 
ber density of H nuclei, and N = 10 22 A^22Cm -2 is the absorbing column density between the nucleus and the gas being con- 
sidered. The M87 nucleus has L x ~ 10 40 ergs s _1 (Wilson & Yang 2002) so that F x = 8.4 x 10 _3 (r/100pc)- 2 ergs cm" 2 s^ 1 . 
If there is 10 6 Mq of gas in a thin disk with outer radius of 100 pc (similar to the fiducial estimates of TB05 and the 
sensitivities being probed by our SMA observations), then the column density normal to the disk would be N22 — 0.3. If 
the vertical scale-height is 10 pc, i.e. for a thin disk aspect ratio of 0.1, then the mean density in the disk is n — 100 cm~ 3 
and the column density through 100 pc of the disk is A22 = 3. For a cloud that is overdense by a factor of 10 from 
the mean, i.e. = 1, we therefore estimate £ c g ~ 1.1 x 10~ 4 . The results of Maloney et al. (1996, their Fig. 3b) 
indicate the gas would be essentially fully molecular, with standard CO abundances, and a temperature of about 20 K. 
If £ e ff were higher by a factor of 10, then the CO abundance is predicted to be reduced by a similar factor. We conclude 
that if the observed 100 pc scale circumnuclear disk in M87 contains a gas mass of ~ 10 6 M Q , allowing gravitational 
instabilities to form over-dense clouds (TB05) , then the present nuclear X-ray fluxes are small enough to allow formation 
of molecular clouds of similar CO content and temperature to giant molecular clouds seen in the Milky Way. As disk 
material approaches the black hole, the increasing X-ray flux will eventually cause the gas to transition to a fully atomic 
state. 

Radiative pumping of lower excitation CO rotational states may also be important in affecting the observed intensity of 
CO(J=2-l) emission (Maloney, Begelman, & Rees 1994). The effectiveness of radiative excitation compared to collisional 
excitation depends on the incident mm flux from the AGN and the gas density. For an observed 230 GHz flux of 1.8 Jy, the 




500 1000 1500 2000 

Velocity (km/s) 

Fig. 4. — Continuum-subtracted spectrum at the positions of the nucleus of M87. The spectrum was extracted toward the continuum peak 
position as shown in Figure 1. The systemic velocity of the M87 black hole is estimated to be at 1290 kms" 1 (Macchetto et al. 1997). The 
regions of the spectrum expected to be line-free are 0-700 kms -1 and 1700-2300 kms -1 , i.e. at least 400 kms -1 from the systemic velocity, 
for which the rms dispersion is 6.8 mjy per 100 kms -1 channel. A few channels in the intervening region show weak emission, but at a level 
too small to be considered a significant detection (see text). 
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flux density at distance r from the nucleus is 4.6 x 10 _13 (r/100 pc) -2 ergs cm" 2 s" 1 Hz" 1 and the corresponding radiative 
excitation rate for J = 1 — » 2 is 7.1 x 10 _8 s _1 . At 20 K the collisional rate coefficient with H2 is 712 = 3.1 x 10 _11 cm -3 with 
a weak temperature dependence. The rates of radiative and collisional excitation will be equal at densities of 2 x 10 3 cm~ 3 . 
Thus molecular gas at densities of n.3 ~ 1 will have its CO excitation temperature raised above the kinetic temperature 
by factors of order unity due to the nuclear radiation source. 

With these caveats in mind, in the following, for simplicity, we proceed using the standard Galactic value to derive H2 
column densities. We also assume that the molecular gas is thermalized up to the J=2-l transition, i.e. that the ratio 
in surface brightness between the 1-0 and 2-1 transition is unity, as is observed in the center of our Galaxy and others 
(Braine & Combes 1992; Sawada et al. 2001). 

3.2.2. Estimates of Molecular Gas Mass 

The continuum-subtracted spectrum in the range 700 to 1700 kms" 1 has a total integrated intensity of 250 Kkms" 1 , 
using a Planck conversion factor of 28.3 K Jy" 1 . Assuming independent la errors of 6.8 mJy per 100 km s _1 channel, 
the lcr error in this total is about 25%. However, as discussed above, the true uncertainties are likely to be correlated 
systematic errors due to continuum subtraction. If the errors are 6.8 mJy per 100 kms" 1 channel and perfectly correlated, 
(i.e. in the same direction for all channels), then the uncertainty in the total intensity is about 80%, i.e. it could be 
consistent with zero. 

Assuming the intensity in the J=l-0 line is the same as that we estimate in the J=2— 1 line, then N(H.2) — 4.5 ± 
(1.1, 3.6) x 10 22 cm -2 , i.e. E = 0.21 ± (0.05, 0.17)g cm" 2 , including He, where the first error assumes independent velocity 
channel errors and the second assumes perfectly correlated velocity channel errors. If this is the mean surface density 
over the 1" wide beam (geometric area = 7r(39pc) 2 ), then the total gas mass in this beam is 4.9 ± (1.2,3.9) x 10 6 M Q . 
Note that because of the disk's inclination to our line of sight, this area would probe regions that are up to 62 pc from 
the black hole. 

If we treat the errors in the 100 kms" 1 channels independently, and integrate over an assumed line width of 500 kms" 1 
(following Braine & Wiklind 1993), then our 3er sensitivity is 2.5 x 10 6 M@, similar to the value they derived, but for a 
region that has an area 170 times smaller and a for a band pass that is almost 4 times wider. 

Another method of estimating an upper limit to the presence of molecular gas within a certain velocity interval comes 
by assuming that three adjoining channels are each required to have a > 3d detection. This would correspond to an 
intensity of 173 Kkms -1 and a mass of 3.3 x 10 6 M Q within a velocity interval of 300 kms" 1 . Requiring > 3cr detections 
in all channels over a 700 km s" 1 velocity interval corresponds to a mass sensitivity of 8 x 10 6 M Q . 

Given the potential presence of correlated systematic errors, our results do not represent a definitive detection of CO 
emission from molecular gas. The data suggest an upper limit on the mass of molecular gas that is present within a 
velocity range ~ 700 km s" 1 about the systemic velocity of ~ 8 x 10 6 Mq. 

To illustrate the nature of a potential emission line of CO, in Fig. [5] we show the expected line profile from a thin 
circumnuclear disk around a 3.4 x 10 9 M Q black hole in M87 with 4.9 x 10 6 Mq of gas uniformly distributed from 50 to 
100 pc and with the disk inclined by 51° to the line of sight and a central velocity of 1290 kms" 1 (Macchetto et al. 1997). 

For a disk fed by Bondi accretion, the size of the disk is limited to within the Bondi radius (about 110 pc, 1.4"for M87). 
This is about the extent of the observed Ha disk in M87 (Ford et al. 1994). The observed velocity gradient of the Ha 
disk is along a position angle approximately perpendicular to the jet axis. However, we do not find evidence for a velocity 
gradient in the potential CO emission features in our data. 

The possibility of seeing CO absorption features against an AGN sub-mm continuum source has been discussed by, 
for example, Barvainis & Antonucci (1994), Maloney, Begelman, & Rees (1994), Wiklind & Combes (1996), and Liszt 
& Lucas (1998). We do not expect local absorption of the nucleus from the 100 pc scale circumnuclear disk because 
of its large inclination angle to our line of sight. Absorption due to CO distributed on larger scales with M87 would 
be expected to produce a relatively narrow (< 300 km s" 1 ) feature centered around the systemic velocity of the galaxy 
(1290 km s" 1 ). This could explain the relative dearth of CO emission that we see from 1000-1300 km s" 1 compared to 
the surrounding channels. Estimating a reduction in flux of ~ 0.015 Jy relative to the continuum level of 1.5 Jy, implies 
r CO(2-i) = 0.01. Assuming this is constant over a 300 km s _1 velocity range and using the conversion factor derived by 
Wiklind & Combes (1996) (for 16 K molecular gas), we estimate a total absorbing column iV(CO) = 1.8 x 10 16 cm" 2 so 
that 7V(H2) ~ 1.4 x 10 20 cm -2 . This is close to, but slightly lower than, the upper limits measured by Braine & Wiklind 
(1993). Thus it is possible that absorption by a larger scale distribution of CO, if present in M87, could have a significant 
influence on the appearance of any CO emission spectrum from the circumnuclear disk, reducing the observed flux near 
the systemic velocity and thus causing its profile to assume a more double-peaked appearance. 

4. CONCLUSIONS 

We have observed 230 GHz continuum emission from the nucleus and jet of M87. The measured fluxes are generally 
consistent with those expected from synchrotron models of these components, although it is possible that thermal dust 
emission makes some sub-dominant contribution in the nucleus. 

The broad spectral coverage of our data allows us to determine the continuum level in regions of the spectrum that 
are expected to be free of line emission from any CO that might be present. After continuum subtraction, we are left 
with a few potential emission features (in 100 kms -1 channels) that are at strengths of about 3-4 times the continuum 
fluctuations in the line-free regions. These features may be due to systematic errors associated with the continuum 
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Fig. 5. — Theoretical (zero noise) spectrum (solid line) of CO(2-l) emission from a thin circumnuclear disk around a 3.4 X 10 Mq black 
hole in M87 with 4.9 X 10 6 Mq of gas uniformly distributed from 50 to 100 pc and with the disk inclined by 51° to the line of sight and a 
central velocity of 1290 kms -1 (Macchetto et al. 1997). The dashed line shows the spectrum averaged into 100 kms -1 wide channels. 

subtraction. Alternatively, if interpreted as CO emission, the implied gas mass is about 5 x 10 6 Mq. Our conservative 
estimate for the upper mass limit of molecular gas that is present in a 700 kms -1 wide velocity range is also at about this 
level, i.e. 8 x 10 6 Mq. Ignoring systematic errors, as is typically done in such estimates, we have a 3a mass sensitivity to 
about 3 x 10 6 M . 

The presence of molecular gas inside the Bondi radius of M87's central black hole is a strong discrimator of different 
accretion models. If present it would indicate that gas has been able to cool to low (T <~ 100 K) temperatures, presumably 
in the confines of a thin accretion disk, whose presence is already indicated by the observations of atomic lines (Macchetto 
et al. 1997). This structure is not part of the paradigm of RIAF models that are fed by Bondi accretion (e.g. Di Matteo 
et al. 2003). Although a RIAF could exist at the center of the disk, much closer to the black hole, the accretion flow at 
its outer boundary would only be that fraction of the mass flux making it through a potentially star-forming disk, not 
the unimpeded Bondi accretion rate. 

The presence of molecular gas is a prediction of the star- forming accretion disk model of TB05, in which strong 
gravitational instability in the disk causes most of the accreted mass to form stars at relatively low-luminosity rather than 
join the black hole. The predicted molecular gas mass is (1 — 5) x 1O 6 M0, which, while quite uncertain, is consistent 
with the observational constraints presented here. Higher sensitivity observations are required to place more stringent 
constraints on the presence or not of molecular gas in the circumnuclear disk of M87. 
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